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ABSTRACT 
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CO 

O \ We calculate the normalized angular three-point correlation function (3PCF), 

q, as well as the normalized angular skewness, S3, assuming the small-angle 
approximation, for a biased mass distribution in flat and open cold-dark-matter 
Oh| (CDM) models with Gaussian initial conditions. The leading-order perturbative 

results incorporate the explicit dependence on the cosmological parameters, the 

^ \ shape of the CDM transfer function, the linear evolution of the power spectrum, 

the form of redshift distribution function, and linear and nonlinear biasing, 

. ^ I which may be evolving. Results are presented for different redshift distributions, 

including that appropriate for the APM Galaxy Survey, as well as for a survey 
with a mean redshift of z ~ 1 (such as the VLA FIRST Survey). Qualitatively, 
many of the results found for S3 and q are similar to those obtained in a related 
treatment of the spatial skewness and 3PCF (Buchalter & Kamionkowski 1999), 
such as a leading-order correction to the standard result for S3 in the case of 
nonlinear bias (as defined for unsmoothed density fields), and the sensitivity 
of the configuration dependence of q to both cosmological and biasing models. 
We show that since angular CFs are sensitive to clustering over a range of 
redshifts, the various evolutionary dependences included in our predictions imply 
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that measurements of q in a deep survey might better discriminate between 
models with different histories, such as evolving vs. non-evolving bias, that can 
have similar spatial CFs at low redshift. Our calculations employ a derived 
equation — valid for open, closed, and flat models — for obtaining the angular 
bispectrum from the spatial bispectrum in the small-angle approximation. 

Subject headings: large scale structure of the universe — cosmology: theory - 
galaxies: clustering — galaxies: statistics 

1. INTRODUCTION 

Characterization of the initial distribution of density perturbations and understanding 
their subsequent evolution into present-day structures are among the central aims of 
cosmology. Particular attention has been focused on the low-redshift distribution of 
matter over large scales, which can probe the linear physics of the early universe and 
therefore test models for the origin of large-scale structure, such as inflation. The n-point 
correlation functions^ (CFs) have become the most widely used statistical tools to quantify 
the observed distribution of matter, in part because they can be easily related to the 
predictions of such models and thus discriminate between them. In principle, recovering 
the current distribution of density perturbations is straightforward; one simply maps the 
three-dimensional distribution of mass with a sufficiently large redshift survey. Existing 
redshift surveys, however, typically contain relatively small numbers of objects, have 
inconvenient shapes, and necessarily suffer from redshift-space distortions, all of which make 
it difficult to obtain precise CF measurements, particularly for the higher-order functions. 

Another possibility is to map the angular distribution of galaxies as projected onto the 
surface of the sky. Angular surveys can sample larger volumes of space and contain greater 
numbers of objects. Furthermore, the effects of redshift distortions need not be considered 
in an angular survey. Thus, if the redshift distribution is well constrained, as is often the 
case for large angular surveys, angular CFs can in practice provide far greater statistical 
power and more insight into the evolution of clustering, than sparser, redshift-survey data 
(Groth & Peebles 1977; Fry & Seldner 1982; Jing, Mo, & Bdrner 1991; Baugh & Efstathiou 



4 The term "correlation function" is understood throughout to refer to the connected, or reduced, 
correlation function. Note that the n-point CFs are valid statistics only if the universe is homogeneous 
on large scales. We invoke the assumptions of large-scale homogeneity and isotropy, so that the CFs depend 
only relative distances and are independent of orientation. 
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1993; Cappi & Maurogordato 1995; Cress et al. 1996; Maddox, Efstathiou & Sutherland 
1996; Gaztanaga & Baugh 1998; Cress & Kamionkowski 1998). 

To date, angular-clustering studies have focussed primarily on the angular two-point 
correlation function (2PCF), zu(9), but with advent of new, deep surveys at various 
wavelengths (radio, optical, IR, x-ray) and more powerful computing resources, increasing 
attention is being focused on higher-order CFs. These are particularly important in 
assessing primordial Gaussianity (Fry 1984; Goroff et al. 1986; Bernardeau 1994; Fry 
& Scherrer 1994; Gaztanaga & Bernardeau 1998; White 1998; Scoccimarro et al. 1998; 
Scoccimarro, Couchman, & Frieman 1998), since to leading order they, unlike the 2PCF, 
reflect the inherently nonlinear gravitational amplification of the initial fluctuations. 
Furthermore, they contain information about the cosmological parameters (Bouchet et 
al. 1992; Juszkiewicz, Bouchet, & Colombi 1993; Gaztanaga, Croft, & Dalton 1995; Bouchet 
et al. 1995; Catelan et al. 1995; Martel 1995; Jing & Bonier 1997; Szapudi et al. 1998; 
Kamionkowski & Buchalter 1998) and about the nonlinear bias between the observed and 
underlying density distributions (Jensen & Szalay 1986; Szalay 1988; Fry & Gaztanaga 
1993; Frieman & Gaztanaga 1994; Fry 1994; Jing 1997; Matarrese, Verde, & Heavens 1997), 
which cannot be obtained from zu(9). 

In particular, the angular three-point correlation function (3PCF), Z(9u,9 23 , 9 3 i), and 
relatedly, the normalized angular skewness, £3(6) (proportional to Z evaluated at zero 
lag, for a distribution smoothed over angular scale O), are the lowest-order intrinsically 
nonlinear statistics and can therefore yield important constraints on models of structure 
formation. Moreover, if different populations are differently biased relative to the underlying 
mass distribution then measuring the angular 3PCFs of these populations can provide 
multiple, complementary constraints. Measurements of s 3 , as well as the normalized angular 
3PCF, q(9 12 , 9 23 ,9 31 ) = Z(9 u ,9 23 , 9 31 )/[w(9 12 )w(9 23 )+w(9 12 )w(9 31 )+w(9 23 )w(9 31 )i 
over quasilinear (QL) scales have yielded results which appear to agree loosely with the 
predictions of the so-called hierarchical model (Peebles 1975; Peebles 1980) for gravitational 
evolution of Gaussian initial conditions (Peebles & Groth 1975; Fry & Seldner 1982; Jing 
& Zhang 1989; Toth, Hollosi, & Szalay 1989; Borgani, Jing, & Plionis 1992; Gaztanaga 
1994; Bernardeau 1995; Cappi & Maurogordato 1995; Gaztanaga, Croft, & Dalton 1995; 
Gaztanaga & Bernardeau 1998; Jing & Bonier 1998; Magliocchetti et al. 1998). It is 
obvious, however, that improved measurements [from surveys such as the Sloan Digital 
Sky Survey (SDSS; Loveday et al. 1998)] and more sophisticated theoretical modeling of 
angular CFs are required. In particular, there has been little theoretical work to date on 
the detailed behavior of angular CFs beyond the two-point function, despite the practical 
advantages they can offer. 
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In this paper, we present a calculation of the full angular 3PCF(], Z, as well as S3, for 
an arbitrary, biased, tracer-mass distribution in flat and open cold-dark-matter (CDM) 
models, assuming Gaussian initial conditions. The calculation, based on leading-order PT 
results, is restricted to the QL regime, where w <C 1, and further invokes the small-angle 
approximation. We take into account such factors as the explicit dependence on the 
cosmological parameters, the shape and linear evolution of the CDM power spectrum, 
and the form of redshift distribution function. We consider both the redshift distribution 
of the APM Galaxy Survey, as well as that of a survey with a mean redshift of ~z ~ 1. 
We also examine the effect of linear and nonlinear bias, which, through an extension 
of the Fry (1996) bias-evolution model to the case of an arbitrary expansion history, 
are allowed to evolve in time. Following BKJ, we define the bias parameters, 6j, by an 
expansion of the unsmoothed tracer-mass field in terms of the unsmoothed dark-matter 
field, rather than by relating the smoothed fields, as is often done (e.g., Fry & Gaztanaga 
1993). While many of the results we obtain are similar to those found in a related paper 
(Buchalter & Kamionkowski 1999, hereafter BK99) dealing with the spatial 3PCF and 
skewness, understanding the detailed behavior of angular statistics in practice can be 
equally, if not more, illuminating. We extend previous derivations of S3 (Bernardeau 1995) 
to include a scale-dependent, leading-order correction which arises in the case of nonlinear 
bias, as defined for the unsmoothed fields, and which becomes large for positive effective 
power-spectrum indices. This behavior, in principle, may allow better constraints on the 
linear- and nonlinear-bias parameters on the basis of large-scale skewness measurements 
alone, or at least differentiate between the smoothed and unsmoothed biasing scenarios. 
For plausible models, we find that s 3 is relatively insensitive to the adopted cosmology, as 
compared with its stronger dependence on the biasing scheme. In general, the presence of 
linear bias tends to flatten and reduce the scale dependence of S3, while a nonlinear bias 
tends to produce a relative increase. The full angular 3PCF, Z, shows similar sensitivities 
to the bias parameters, but can also depend significantly on the cosmological model, and in 
particular on the shape of the power spectrum. We find that the configuration dependence 
of Z is in general more complex than the simple hierarchical model suggests, and, since 
angular statistics probe clustering over a range of redshifts, is also sensitive to evolutionary 
effects. These properties can be used to discriminate between various models which can 
yield nearly degenerate predictions for s 3 , or for the spatial 3PCF at z — 0, such as evolving 
versus non-evolving bias models, or open versus flat cosmological models. Neglecting the 
full geometric and evolutionary variation of Z throws away valuable information which 
might, in practice, be harder to obtain from three-dimensional data. In several instances, 



5 Frieman & Gaztanaga (1999), using a different method, have independently derived results for Z. We 
have compared our results with theirs and find that they agree for some of the models we have considered. 
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we illustrate how accounting for time dependences in, e.g., the power spectrum or the bias, 
can impact the resulting calculations, especially in the case of high-redshift surveys. Our 
theoretical results, in general, should be comparable to current and future measurements of 
angular clustering. 

In §2 we outline the calculations used and present results for the normalized angular 
skewness. In §3 we derive predictions for the full angular 3PCF and in §4 we summarize 
our conclusions. Since the calculations are greatly simplified by working initially in Fourier 
space, we first obtain the angular bispectrum — the two-dimensional projection of the spatial 
bispectrum — and then derive the angular 3PCF from this. A derivation of the equation 
that relates the spatial and angular bispectra is presented in the Appendix. This result is 
essentially a Fourier-space generalization of Limber's equation for the 3PCF. 



2. ANGULAR SKEWNESS 

Since angular CFs involve projections of their full spatial counterparts, our derivations 
must involve integrals over the line of sight, which will in general depend on the adopted 
cosmological model. We shall consider open and flat Friedmann-Robertson- Walker (FRW) 
cosmologies with a possible cosmological constant, so that the scale factor, a(t), satisfies 
the Friedmann equations, 

- = H E(z) = H Jq (1 + zf + (1 - tt - n A ){l + zf + fi A , 
a v 

- = H 2 [Q A -n (l + z)y2}, (1) 
a 

where Qq is the present nonrelativistic matter density in units of the critical density, Q\ 
is the contribution of the cosmological constant to the total present energy density, H is 
the present value of the Hubble parameter, and a dot denotes a derivative with respect to 
time. The scale factor is chosen such that aoH = 2, where H = 100 h km s^ 1 Mpc -1 . If 
we take our position as the origin, w = 0, then the angular-diameter distance to an object 
at redshift z is given by 

where S(x) = {sinhx, x, sin a;}, respectively, for open, flat, and closed geometries. For 
an Einstein-de Sitter universe, w(z) = f(z). The distance to the horizon is given by 
i] = w(oo). 
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We wish ultimately to derive an expression for the angular 3PCF of a 
distribution of tracer masses (e.g., galaxies, quasars, clusters, radio sources, etc.), 
Z (#12, #23, 6*3i ) = (p(6 1 )p(O 2 )p( y 03)) , where represents two-dimensional coordinates on the 
sky, p(0) = — is the fractional perturbation to the unsmoothed tracer-mass 

density field, S(0), and angular brackets denote an average over direction in the sky, with 
fixed values for the distances % —\0i — 0j\. Though statistical homogeneity and isotropoy 
ensure that the 3PCF will only depend on the relative angular separations between the three 
points, measuring Z for all possible geometric configurations can nonetheless be a daunting 
task. For this reason, many studies have focussed instead on the normalized angular 
skewness, s 3 = (pq(0)) / {p%>{0)) , obtained from the moments of counts in cells smoothed 
with an effective angular scale half-width (it is implicitly assumed that £ oc n, where n is 
the number of discrete counts). This one-point statistic is more readily calculable than the 
full angular 3PCF, while still preserving information about the overall scale-dependence of 
Z. 

In deriving the following results, we implicitly make use of the small-angle 
approximation in assuming that a small patch of sky can be treated using a Fourier 
expansion (Peebles 1980). In this approximation, the component of the three-dimensional 
wave vector along the line of sight is taken to be negligible compared to the orthogonal 
components, so that the angular CFs depend only on the latter. We later comment on this 
assumption. We define pe using a two-dimensional spherical top hat window function, so 
that in terms of the Fourier components, p(n), of the unsmoothed two-dimensional density 
field, we have 

Pe{6) = I ^^K^M^e), (3) 

where k is a two-dimensional wave vector with magnitude k, W(x) = 2Ji(x)/x is the 
Fourier transform of the two-dimensional spherical top-hat window function, and is an 
angular smoothing radius. We can then obtain the angular 2PCF, as measured using counts 
in cells, 

t*e(0i 2 ) = <Pe(*i)pe(*2)> = / / ^y^ e^^^W^W^e) (p>i)p> 2 )) ■ 

(4) 

We define 

(p(ki)p(k 2 )) = (2n) 2 5 D ( Kl + K 2 )P P («), (5) 

where the angular, or projected, power spectrum of the tracer mass, P p (usually written 
in terms of C/s in papers on the cosmic microwave background), can be related to the 
full tracer-mass power spectrum, P(k,w), by a Fourier-space analog of Limber's equation 
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(Kaiser 1992), 

Here, dN/dw is a w-space selection function for the tracer mass, normalized such that 
Jo* dw{dN/dw) = 1. 

We assume the fractional perturbation to the unsmoothed three-dimensional tracer 
mass density field, 6(r), may be expanded in terms of the local perturbation to the 
unsmoothed, underlying matter field, 6 m (r), via 

5 = b l 5 m + b j6 2 m + ---, (7) 

where h\ is the linear bias term, 62 the first nonlinear term, etc., and 5 m is itself written as 
5 m = 5$ + 6$ + ■ ■ ■ , where 6j%' <C $m~ l \ $m ls the linear solution, and 5$ characterizes 
the leading-order departure from the Gaussian initial conditions. The linear solution for 
the underlying spatial density contrast has the separable form 

5^(r,w) = D(w)5^(r,0), (8) 

so that fluctuations evolve simply as the linear growth factor, D(w). This in turn implies 
that, to leading order, the spatial and time dependence of the power spectrum can likewise 
be factorized. For FRW cosmologies, the growth factor is given (as a function of redshift) 
by 

For an Einstein de-Sitter universe D(w) is simply the scale factor, a(w). Assuming the 
unsmoothed three-dimensional tracer-mass density contrast to be related to the underlying 
distribution via equation (|7]), we can write the leading-order result for the full linear power 
spectrum of the tracer mass, 

P(k,w) = AblD 2 {w)k n T 2 {k), (10) 

where k n is the primordial power spectrum, A is the overall amplitude, and T(k) is a 
model-dependent transfer function. Substituting equations (|TUD, (||), and (^) into (|J), 
taking Oi = 6 2 = 6, and x = kQ, we arrive at the angular variance (the 2PCF at zero lag), 

«*«» = (pim) = ^A- £ dw 2 b .Lim jT dx x ^w\ x) T^ie w) , (id 

where the 61 term is left inside the w integral to allow for the possibility of bias evolution. 
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The counts-in-cells angular 3PCF is given by 

x (p(Ka)p(K 2 )p(K 3 )) vv(Kie)w(K 2 e)w(K 3 e). (12) 

The angular, or projected, bispectrum, is denned via 

(p(k 1 )p(k 2 )p{k 3 )) = (2tt) 2 5 d (ki + n 2 + k 3 )S p (ki, k 2 , K3), (13) 
and is related to the full three-dimensional bispectrum, B, by 

cfay fdN\ 3 

B p {k\, «2, « 3 ) = y o ^ll"^") B(k 1 /w,k 2 /w,k 3 /w,w) (14) 
(see Appendix). The full bispectrum is given, to leading order in PT, by 



B(ki,k 2 ,k 3 ,w) = P(kx,w)P(k 2 ,w)l — 



(k\ k 2 \ 2 
1 + ix + cos-0 - — h — + (1 — n) cos V 

\K2 Kl/ 



+ ^| + (eye), (15) 

(Fry 1984; Goroff et al. 1986; Matarrese, Verde, & Heavens 1997) where if) is the angle 
between k\ and &2, and \l is a function of the expansion history, equal to 3/7 for an 
Einstein-de Sitter universe (Peebles 1980) and differing from this value only slightly for 
other reasonable choices of the density parameters (Bouchet et al. 1992; Bernardeau 
1994; Bouchet et al. 1995; Catelan et al. 1995; Martel 1995; Scoccimarro et al. 1998; 
Kamionkowski & Buchalter 1998). Note that the expression for B includes the dependence 
on b 2 , since the leading-order result for the 3PCF includes second-order terms in the 
expansion of equation (0). 



Evaluating equation ( |i~2"D with Q\ = 6 2 = 6 3 = 6 } X{ = kB, taking m to lie in the 
if) — direction, and substituting equations (0), QI4]) , and (|j"5|) we obtain, after a little 
algebra, an expression for the area-averaged skewness, 

3A " r A (dN\ 3 b^D\w) 



(27r) 3 e 2 "+ 4 Jo \dw J w 2n + 4 

x 



rOC roo 

/ dx 1 x , l +1 W(x 1 )T 2 (x 1 /ew) / dx 2 x^ +1 W(x 2 )T 2 (x 2 /Qw) 
Jo Jo 

1 

h 



x r d<f> W(x 3 ) ( 1 [(1 + //) + cos ^(- + ^) + (1 - //) cos 2 J + ps) (16) 
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where one integral vanishes under the requirement n 3 = + k 2 ), and a factor of 3 arises 
from symmetry considerations applied to the 2 cyclic permutations in equation (|15|). Noting 
that £3 = \jx\ + x\ + 1x\x 2 cos ip, we can evaluate the ip integrals by using the summation 
theorems for Bessel functions (Gradshteyn & Ryzhik 1980; Bernardeau 1995) which yield: 

sin 2 ^ W(x 3 ) = tcW{x 1 )W{x 2 ), (17) 







r Zn x<> 

/ di/> (1 + — cost/0 W(x 3 ) = 2tcJ {x 2 )W{x 1 ), (18) 

JO X\ 

/ dipW(x 3 ) = 27r^(2j + l)>V 2 , + i(a; 1 )W 2 , +1 (a; 2 ), (19) 
Jo 3=0 

where W n (x) = 2J n (x)/x, so that W(x) = Wi(x). Using the above results gives 



3A 2 no ( dN\ 3 bi i D i { 



Pe(0)) = T^E^Ti I dw — 



w 



(27r) 2 e 2 ™+ 4 Jo \dw J w 



2n+4 



/ dx 1 x1 +L W(x 1 )T 2 (x 1 /Qw) / dx 2 x 2 i+i W(x 2 )T 2 (x 2 /<dw) 
Jo Jo 

|^ [Jo(x 2 )Mx 1 ) + J (x l )W(x 2 )} - l^Ay V (x 1 )W(x 2 ) 

b 00 1 
+3 E(2J + l)W 2j+1 (x 1 )W 2 , +1 (x 2 ) . (20) 
°i j=0 J 

We re-emphasize that our results in general are based on leading-order PT and are 
thus restricted to the QL regime, i.e., scales large enough so that the rms density contrast 
fluctuations are small compared with unity. Numerical simulations and observations both 
confirm that higher-order nonlinear corrections have little impact on spatial statistics on 
QL scales (Szalay 1988; Toth, Hollosi, & Szalay 1989; Gott, Gao, & Park 1991; Fry, Melott, 
& Shandarin 1993; Jain & Bertschinger 1994; Scoccimarro et al. 1998). By comparing 
with iV-body results, Gaztanaga & Bernardeau (1998) test the validity of PT as applied to 
angular statistics and find good agreement on scales 6 > 1°, with the details depending on 
the shape of the power spectrum. In addition, simulations show that QL scales can still 
obey leading-order PT even when smaller scales have become fully nonlinear (Bouchet & 
Hernquist 1992), and further that the predictions of QL PT may hold even on scales where 
the rms fluctuation exceeds unity (Bernardeau 1994; Baugh, Gaztanaga, & Efstathiou 1995; 
Fry, Melott, & Shandarin 1995; Szapudi et al. 1998). 

Our results have also relied on the use of the small-angle approximation. The Fourier 
integrals used to calculate w and Z, however, extend over all values of n, including 
small values for which this approximation breaks down, but which nonetheless may 
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contribute significant weights to the integrands on degree scales.^ Moreover, the small-angle 
approximation might appear to go against the assumption of quasilinearity; in practice, the 
degree to which these two assumptions are mutually plausible will in general depend on 
the selection function in question. Despite these considerations, our results for S3 and q on 
scales near 1° do agree well with data from observations (e.g., Gaztanaga 1994; see Figure 
0) and iV-body simulations (Frieman & Gaztanaga 1999) on these scales. Bernardeau 
(1995) investigates the validity of the small-angle approximation for models with different 
redshift selection functions, dN/dz, and finds good agreement with numerical results for 
S3 on scales near and below 1°. Gaztanaga & Bernardeau (1998) find that the small-angle 
approximation can yield reasonable agreement with numerical results for S3 out to the 
5° — 10° scale. Verde et al. (1999) perform an exact calculation for the angular 3PCF using 
a spherical-harmonic decomposition, and a more complete comparison and assessment of 
the small-angle results against the exact results are presented therein. 



2.1. Effects of the Power Spectrum 



Equations (|TT| ) and (|20| ) are one-point results which, for a given choice of fi > ^A> b±, 
and 6 2 , depend only on the smoothing radius, 0, and the form of the power spectrum. If we 
assume, for computational simplicity, a scale-free power spectrum, P{k) oc k n [T(k) = 1], 
then we can calculate an expression for S3 where the x integrals can be separated from the 
w integration and evaluated analytically; combining equations ( |i~l~l) and (|20|) , we obtain 



s 3 



where 



and 



(g|(g)> 

(pUo)) 2 




36 9 (1 

1 \-u-l 

7 14 V3 P 



3^n 
6 b\ [i 



A(n)] 



R 3 (n) 



JS° dw (dN/dwf w-^ 2n+ ^D\w) 



JS° dw (dN/dw) 2 w-( n+2 W 2 (w) 



2 ■ 



A(n) 



E£i(2j + 1) U^dxx^J^J^x)]' 



(21) 



(22) 



(23) 



[^dxx^Jl{x)f 

Note that while the term in curly brackets in equation ( |2"T| ) depends only very weakly 
on fl through fi, the i?3(n) term depends on the cosmological model both through the 
presence of the linear growth factor, D (which is often neglected), as well as through the 
explicit dependence of w on the expansion history [see equation @]. Equations (^TT)-(|2 



6 We note that this difficulty could be avoided if one considers the bispectrum directly. 
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Fig. 1. — The variation of the A term in equation with spectral index n. Though 
negligible for n in the range — 2 < n < 0, this term rapidly becomes significant for n > 0, 
due to the impact of small-scale fluctuations. 
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are valid for n in the range — 2 > n > 1; the individual terms in the series for A(n) can be 
evaluated explicitly (Watson 1966) and we find that A(n) <C 1 for —2 < n < 0, but diverges 
for n > 0, as illustrated in Figure [I]. 

The bracketed n-dependent terms in equation (|2l|) arise from the effects of smoothing 
the contiguous field, 5(r), defined in equation (|7]). In particular, since smoothing at a fixed 
physical mixes different mass scales, the higher-order moments for smoothed fields will 
generally contain additional terms reflecting the scale dependence of the rms fluctuation 
(Bernardeau 1994, 1995). The n = —2 case thus recovers the no-smoothing result for 
S3 (Bernardeau 1995), since the rms fluctuation in this case is independent of scale [see 
equation ([□])]. Though the filter is designed to separate out the nonlinear, small-scale 
fluctuations, an unbounded scale-free power spectrum with n > produces so much 
small-scale power that the fluctuations induced by smoothing over nonlinear scales become 
comparable to or greater than the large-scale linear perturbations, leading to the divergence 
in A(n) seen in the Figure. 

Equation ( pip differs from the known result for S3 assuming constant, nonlinear bias 
(Bernardeau 1995) only in that it contains the A(n) term. The difference is attributable 
to the fact that previous authors usually define the bias parameters using the form of 
equation (0), but relating the smoothed, rather than the unsmoothed, observed and 
underlying density fields. Such a definition is technically nonlocal, only fixing the linear- 
and nonlinear-bias parameters at the chosen smoothing scale. While it can be shown that 
the "smoothed" and "unsmoothed" linear-bias parameters are identical, the value of the 
smoothed nonlinear-bias parameter will, unlike the unsmoothed 62 defined in equation (^) 
in the limit of continuous fields, depend on the smoothing scale. Comparing the result for 
s 3 with those of Fry & Gaztanaga (1993) and Bernardeau (1995), we can infer that the 
smoothed nonlinear-bias parameter will vary as ^[l + A(n)]. In the limit of no smoothing 
(n = —2), the two parameters are identical, as expected, but can differ dramatically on large 
scales (see §2.2). Similar correction terms, arising from nonlinear bias, are expected for 
higher-order angular moments as well. Such terms might be used to discriminate between 
the smoothed and unsmoothed bias pictures, and possibly distinguish between linear and 
nonlinear bias based on the skewness alone, as discussed in §2.2 and in BK99. Unless stated 
otherwise, it will be hereafter understood that "bias" shall refer to the unsmoothed bias 
prescription defined by equation (|7|). 

While the power spectrum for any viable cosmological model will not be given by a 
simple power law, it can be shown, using the results of Bernardeau (1995) and Gaztanaga & 
Bernardeau (1998), that an exact result for S3 can be obtained for arbitrary power spectra, 
using a properly redefined effective index (Scoccimarro 1998). This result, analogous to 
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that in the spatial case (Bernardeau 1994), is strictly valid, however, only for linear biasing, 
since it fails to account for the above-mentioned large-scale variation associated with 6 2 - 
Since we wish to consider nonlinear biasing, and because, unlike S3, no corresponding exact 
result can be obtained for the full angular 3PCF in the case of arbitrary power spectra, 
even for linear biasing, the results of this paper are based entirely on numerical integration 
of the appropriate equations. Many of these will involve highly oscillatory integrands, 
making them difficult to evaluate; we have performed checks on the numerical accuracy of 
our results, and found them to be good to within a few percent. 



For realistic CDM models, we employ a transfer function given by 

m(l + 2.34g)/(2.34g) 
[1 + 3.89g + (16. lg) 2 + (5.46g) 3 + (6.71g) 



r<«) = :::::: Try. „ l/t m 



(Bardeen et al. 1986), where q = k p /Y, V « Q Q h, and k p is the physical wavenumber in units 
of h Mpc -1 , related to the comoving wavenumber, k, through our adopted normalizations 
by k p = fc/6000 h Mpc -1 . The result for s 3 is then obtained by using equation ( f2"4"|) together 
with ([Tip and fl2"0"|). Note that in the case of a scale-dependent transfer function the x and 
w integrals cannot be separated. 

Our calculations also require the assumption of a selection function (SF) along the line 
of sight, and we consider two functional forms chosen to characterize low-redshift (z ~ 0.1) 
and high-redshift surveys (z <~ 1). For the former, we consider the SF of the Automatic 
Plate Measuring (APM) Galaxy Survey, 



dN 2 
— — oc z exp 

dz 



1.412z c , (25) 



where z m is the median reshift of the survey, taken to be 0.12. This functional form 
provides very good fit to the APM redshift distribution over the entire APM magnitude 
range (Baugh & Efstathiou 1993; Maddox, Efstathiou, & Sutherland 1996; Gaztanaga & 
Baugh 1998), and should also approximate comparable optical surveys, such as the SDSS. 
To compare the expected results for a survey such as the APM, which reaches only modest 
redshifts (z < 0.4), with those from a much deeper angular survey, we also employ the 
w-space SF of Kaiser (1992) (intended to mimic a magnitude-limited survey), 



dw wl+ a T[(l +a)/0\ 



w, = w(T[(2 + a)/P}T[(l + a)/P}), (26) 



where w is the mean conformal lookback time. We will hereafter choose a = 4, (3 = 4, and 
w = 0.35, which yields a peak redshift of z ~ 0.8, so that this model, which we will refer 
to as the "high-;?" SF, might characterize, for example, a radio survey such as the VLA 
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Fig. 2. — Redshift-space plot of the APM (z m = 0.12) and high-z (z ~ 1) selection functions 
employed in the calculations. For comparison, we also display the APM SF assuming 
z m = 1.0; both this and the high- z SF have been assigned three times the normalization 
of the low-redshift APM SF for the purposes of representation. 
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FIRST Survey (Becker, White, & Helfand 1995), with a median redshift of z m ~ 1 (Cress 
& Kamionkowski 1998). Figure ^| shows the above APM SF together with the high- z SF, as 
well as the APM SF obtained using z m = 1.0, for comparison. Baugh & Efstathiou (1993) 
and Maddox, Efstathiou & Sutherland (1996) find that angular statistics in the APM 
Survey depend much more sensitively on the peak redshift of the SF, rather than on its 
precise shape; we find this conclusion to hold in general for SFs we consider. Bernardeau 
(1995) finds that the small-angle approximation works reasonably well for an APM-like SF; 
for deeper surveys in which a fixed angle can probe larger scales, this approximation is 
expected to provide still better results. 

In Figure ^, we plot the variation of S3 with B, assuming the APM SF, for an unbiased 
{pi = 1, 6 2 = 0) tracer- mass population in four flat (fl + Q\ = 1) cosmological models: 
standard CDM (SCDM; Q = 1, h = 0.5, n = 1), a tilted CDM model (TCDM; Q = 1, 
h = 0.5, n = 0.8), CDM with a cosmological constant (ACDM; Q = 0.4, h = 0.5, n = 1), 
and CDM with a high Hubble parameter (HCDM; Qq = 1, h = 0.75, n — 1), as well as 
an open CDM model (OCDM; Q = 0.4, h = 0.5, n = 1). The thin horizontal line in the 
Figure shows the (virtually f2 -independent) semi-analytic results obtained for a scale-free 
power spectrum with a canonical index value of n = —1.2 (Peebles 1980). The results for S3 
in general are very similar to those found by BK99 for the normalized spatial skewness. The 
scale-dependent CDM transfer function naturally induces a dependence on the smoothing 
scale, 6, and compared to equation (|2"T|), a more substantial dependence on the combination 
of cosmological parameters T = Q h. Yet, although the five cosmological models shown in 
the Figure span a broad range of parameter values, 0.2 < T < 0.75, they yield fairly similar 
results, and are certainly all consistent within the limits of current observations (see Figure 
|7|). The slight variation seen is mostly due to the dependence on T; there is little difference 
between the SCDM and TCDM models, which differ only in the value of n, and virtually 
no difference between the OCDM and ACDM models. If the analysis is restricted to the 
range 0.2 < T < 0.3, as suggested by current data (Bartlett et al. 1998), one can infer from 
the Figure that the predicted values for S3 for a survey such as the APM are relatively 
insensitive to the adopted cosmological model. Figure f| shows the results for 53(G) for the 
same five cosmological models, assuming the high- 2 SF. The primary dependence is again 
seen to be with T, and while there is a somewhat greater spread between the different 
models, the differences between realistically viable models are still relatively small compared 
to those induced by the bias dependences which we explore in the next section. 

Overall, we find that the high- z SF curves for these models show greater variation 
of S3 with O than their APM SF counterparts. This is simply due to the fact that in a 
deeper angular survey, a fixed angle probes clustering over a broader range of physical 
scales, including very large scales, while conversely, in a shallower survey, a fixed angle 
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Fig. 3. — The variation of S3 with 0, for the APM SF, in the five cosmological models 
described in the text, assuming no bias. Note that the differences between these models 
are relatively small. The large values of S3 at high are due to the shallowness of the SF. 
The solid, horizontal line shows the result obtained from equation (|2l|) for a scale-free power 
spectrum with n = —1.2. 
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Fig. 4. — Same as Figure ^ but using the high-z SF. The models differ slightly more than 
in the APM case, and attain smaller values of S3 at high since these angular scales, in a 
deeper survey, incorporate larger physical scales where the effective power-spectrum index 
is larger. 
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is sensitive to clustering only over a relatively narrow range of smaller physical scales. 
Thus, at larger values of 9, the high-z result incorporates physically large scales and tends 
toward smaller values for S3. These can be roughly approximated using equation (^Tj) 
by s 3 « R 3 (n k ) x [36/7 - (3/2) (n* + 2)] « R 3 (n k ) x (9/14), where n k , the value of the 
effective spectral index of the (untilted) linear-theory CDM power spectrum, tends towards 
unity at large scales (recall that R3 is of order unity and independent of 0). For the APM 
case, however, S3 remains higher at larger 0, since even large angular scales in a shallower 
survey still correspond only to smaller physical scales, where the effective power-spectrum 
index is negative. At smaller 0, both sets of curves tend toward the expected larger values 
of S3 ~ Rz{n k ) x (93/14) obtained with n k m —3 — the effective power-spectrum index at 
small scales; for tilted models, S3 tends toward a value n k — 1 smaller. Note that in both 
Figures, the semi-analytic results for a scale-free power spectrum with a canonical index 
value of n = —1.2 (thin, horizontal lines) provide poor fits to the CDM predictions. Other 
authors have investigated QL results for the normalized angular skewness using CDM 
models (Gaztanaga & Frieman 1994; Frieman & Gaztanaga 1994; Gaztanaga & Baugh 
1998; Gaztanaga & Bernardeau 1998), and find similar results. 



2.2. The Effects of Bias and its Evolution 

The results of the previous section already allow for the existence of a constant, 
nonlinear bias between the tracer-mass and the underlying matter distributions, but can be 
generalized by allowing the bias parameters to evolve with time, as is suggested both by 
theory and observations (Fry 1996; Peacock 1997; Matarrese et at 1997; Steidel et at 1998; 
Cress & Kamionkowski 1998; Catelan et al. 1998; Taruya, Koyama, & Soda 1998; Tegmark 
& Peebles 1998; Colin et al. 1998; Baugh et al. 1998; Magliocchetti & Maddox 1998). 
Fry (1996) proposes a bias-evolution model, which assumes that objects in an Einstein-de 
Sitter universe form at a fixed formation redshift, Zf, by some arbitrary local process which 
induces a bias at that epoch, and are subsequently governed purely by gravity. We note 
that this model fails to account for merging, and cannot produce an anti-bias (pi < 1). 
BK99 discuss a generalization of the Fry model to the case of arbitrary expansion history. 
Other, more general models have been proposed (e.g., Tegmark & Peebles 1998), but we 
employ the generalized Fry model as a first approximation to investigate the effects of bias 
evolution, and merely quote the relevant results below. The bispectrum in this model is 
given by 



C x (w) + C 2 H cos V> ( ^ + j + C 3 (10) cos 2 v 
\k 2 fa) 



B(fa,k 2 ,k 3 ,w) = P(fa,w)P(k 2 ,w) 

+ (eye), (27) 
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where 

Ci{a) = wi m ^K-mw-2 m+ K t (28) 

" nM-r (29) 

+ fax. - If 

rf(tu) = D(w)/D(w*), and a subscripted asterisk denotes the value of that parameter at 
the epoch of formation. Note that we have now ignored the very weak dependence of the 
bispectrum on /x; the dependence on the expansion history, i.e., on the species contributing 
to the total energy density, is contained in the growth factor, D(w). The result for S3 can 
then be extended to the case of evolving bias by using the bispectrum of equation fl27|) in 
equation (|16|), which amounts to making the substitutions 

I— >C7 8 ( W ), ^ — > dW - 2C 2 (w) + C 3 (w) (31) 

Ol Oi 

inside the brackets in equation (|20|), and replacing b\ by 1/C 2 (w). 

For the APM and high- z SFs, respectively, Figures [| and ^ display the results for 
53(0) in five different bias scenarios, each employing the SCDM, ACDM, and OCDM 
models above: an unbiased scenario (61 = 1, &2/&1 — 0, taken from Figures |3] and |4j), 
a non-evolving, linear bias (61 = 2, b^jb\ = 0), a non-evolving, nonlinear bias (61 = 2, 
&2/&1 — 1), an evolving, linear bias = 5, 62* /^i* = 0); an< ^ an evolving, nonlinear bias 
(&i» — 5, &2*/&i„ — 1)- For the non-evolving bias scenarios, the ACDM and OCDM models 
are taken to be identical, as suggested by Figures ^ and f|, and therefore only the ACDM 
models are plotted in these cases. These two models do, however, differ slightly more in 
evolving bias scenarios, due to the different time dependences of their linear growth factors, 
and the differences between them would increase with decreasing Qq. For the evolving cases 
we hereafter assume a formation redshift of Zf = 5. 

The qualitative dependences of S3 on the bias scenario are similar to those of its spatial 
counterpart (BK99). Specifically, it is clear from these Figures that the dependence of the 
angular skewness on the biasing scheme can potentially be far more significant than that 
upon the cosmological parameters within a given scheme. For example, adding a very small 
linear bias to the low-r models would yield a result very similar to the unbiased SCDM 
predictions. Generally, the presence of any significant linear bias, b± > 1, tends to reduce 
the dependence of S3 on 0, as compared with the unbiased cases (solid lines). Furthermore, 
an observed S3 curve which is far below the predicted unbiased result can only be achieved 
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Fig. 5. — Plot of Ss(Q), assuming the APM SF, for the five bias scenarios described in 
the text, each employing the SCDM model (thick lines), ACDM model (medium lines), and 
OCDM model (thin lines). For the non-evolving bias cases, the latter two models are taken to 
be identical, as suggested by Figure [3|. Note the strong sensitivity to the biasing scheme and 
the presence of the upturn in S3 at large G in the nonlinear bias cases. As expected, linear 
biasing tends to lower the predicted curves, while nonlinear biasing tends to raise them. The 
evolving bias models yield less dramatic shifts than their non-evolving counterparts, since 
evolving bias can effectively act as a constant bias, which at late times will be small. 
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Fig. 6. — Same as Figure |5|, but using the high-z SF. Note again the impact of the biasing 
scheme, the signatures of linear and nonlinear biasing, and the upturn in S3 seen in cases 
where 6 2 7^ 0. the Figure also illustrates one example of the possible degeneracies of s 3 
with respect to the cosmological parameters, the bias parameters, and the evolution in the 
latter. In this case, the predictions for the SCDM and ACDM models with evolving linear- 
bias (fei, = 5; short-dashed lines) are respectively similar to the SCDM and ACDM/OCDM 
models with constant, linear bias (bi =2; dotted lines). 
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by a large linear bias term (assuming 6 2 > 0), while one well above this value can only arise 
from the existence of a non-zero 62 term or from anti-biasing (pi < 1). The curves for an 
evolving bias can produce less drastic shifts than their non-evolving counterparts, despite 
the fact that the respective bias terms are initially larger. This is because the evolution 
towards an unbiased state can effectively mimic a smaller (especially at late times), constant 
bias, as demonstrated by BK99. 

Again we find that the APM curves in Figure [| are flatter than their high-z counterparts 
in Figure |6], for reasons explained in §2.1. Note also that models with nonlinear bias exhibit 
an upturn at large 0. The upturn arises from the contribution of the higher-order terms 
in the sum in equation (|2Up, which, in the presence of a non-zero b<i term, can become 
significant on large scales, where the effective power-spectrum index becomes positive (see 
Figure 1). BK99 perform the analogous three-dimensional calculation and find similar 
large-scale behavior of the normalized spatial skewness, 5*3. Depending on the cosmological 
model, this transition here occurs at angular scales of about a few degrees for the APM 
SF, and slightly higher values in the case of the high-;? SF. Thus, the SCDM model, whose 
power spectrum peaks at larger wavenumber than those of the low-r models, will exhibit the 
upturn at smaller values of G, while the low-r models with nonlinear bias are expected to 
show an upturn only at larger angular scales, where they begin to sample positively-sloped 
regions of the power spectrum. However, unlike the case for S3 (BK99), the exact location 
and severity of this upturn in S3 for the different cosmological models is influenced not only 
by the shape of the power spectrum, but also by the weightings arising from the SF (see 
Figure |2|) and from the terms proportional to powers of w in equations (^) and (|TT|) . For 
example, the APM SF, which resembles a low-to spike, exacerbates the upturn by effectively 
sampling the w integrals only around the low peak value of w(z m ), where the w~^ 2n+4 ^ 
factor is extremely large. In the case of the high-z SF, while large values of still sample 
a positively-sloped region of P(k), the SF itself, which is everywhere much shallower than 
in the APM case, is effectively zero where the w~( 2n+4 ) factor is large, and peaks at a high 
value of w, where this factor is very small. As seen in Figure |6], this behavior tends to lend 
less support to the divergence of the sum in equation fl2"0p at these scales, even in the case 
of a large, constant value of 62- Note that for both SFs, the non-evolving, nonlinear cases 
exhibit more of an upturn, since in these cases the large value of 62 is maintained for all 
time, whereas in the evolving cases, the divergence is less severe due to the fact that the 
tracer population evolves towards an unbiased state with time. 

The upturn in S3 at large O in nonlinear-bias models differs from the decreasing 
behavior expected if the higher-order terms in the sum in equation (|20| ) are neglected, i.e., if 
the bias parameters, as defined for smoothed density fields in equation (^), are are assumed 
to be constant, as in peak-biasing or halo-biasing models (Mo & White 1996). BK99 show in 
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the case of the normalized spatial skewness that this different large-scale behavior might be 
used to obtain better constraints on b\ and b 2 on the basis of skewness measurements alone, 
since combinations of linear and nonlinear bias which yield effectively degenerate results for 
S3 on smaller smoothing scales, can, at least in principle, be distinguished by considering 
smoothing scales where the relevant effective index (Bernardeau 1994) is positive, i.e., 
scales greater than 100 h^ 1 Mpc in typical CDM models (see Figures 2 and 5 in BK99 and 
related discussions). Similarly, we can infer from Figures [5] and |]that combinations of these 
parameters which are degenerate at small G become more distinct at large (due to the 
upturn), and that comparing QL PT predictions with large-scale measurements (0 > 5°) 
of S3 might also resolve this degeneracy. If this large-scale behavior is not seen, this would 
argue in favor of smoothed-bias prescriptions, rather than the unsmoothed formulation we 
employ. At such scales, however, one might be limited by the variance induced from the 
smaller number of independently sampled cells (Seto & Yokoyama 1998), and may require a 
more accurate derivation which does not assume the small-angle approximation, but instead 
employs a full spherical harmonic decomposition (Verde et al. 1999). 

Figures |5] and || also provide illustrations of one of the potential degeneracies that can 
arise when considering the normalized skewness. It is known that the one-point statistics S3 
and S3 can be nearly degenerate with respect to combinations of cosmological and biasing 
models, with respect to constant b\ and 62 (barring large-scale measurements, as mentioned 
above), and also with respect to evolving vs. non-evolving bias models (BK99). This latter 
degeneracy is illustrated in Figure |5] by the almost identical predictions of the OCDM model 
with evolving, linear bias and the OCDM/ACDM model with constant, linear bias, and in 
Figure || by the similarity between the SCDM models with constant, linear bias and with 
evolving, linear bias. It would clearly be impossible to distinguish between such similar 
predictions given present data of the type shown in Figure [F[ 

As an example of what can be learned from current measurements, we compare in 
Figure [7| our predictions for S3 using the APM SF with data from the APM Galaxy Survey, 
containing over 1.3 x 10 6 galaxies (see Gaztanaga 1994 and references therein). The data 
points (Gaztanaga 1997) are not corrected for the effects of source fragmentation, and the 
error bars shown are derived from the scatter in four sub-zones, plus shot-noise corrections, 
and are thus extremely conservative, especially at large 0. The curves plotted are for 
the SCDM and ACDM models, with "best-fit" values for 61 of 1.14 and 1.22, respectively 
(assuming only constant, linear biasing). It should be noted that discrepancies between 
the data and our QL PT predictions are expected at small scales, due to: 1) shot-noise 
fluctuations arising from the discreteness of the observed distribution [i.e., the (p'oj are no 
longer reliable estimators of the area-averaged CFs (Peebles 1980; Gaztanaga 1994, 1997)], 
2) estimation biases arising from the fact that the ratio of two unbiased estimators is not 
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Fig. 7. — Comparison of APM-SF predictions for S3 in the SCDM and ACDM models 
(with linear bias), with data from the APM Galaxy Survey (Gaztanaga 1994, 1997). While 
deviations are expected at small scales (from nonlinear and discreteness effects) and at large 
scales (from sample variance and the breakdown of the small-angle approximation), the data 
agree reasonably well, in the QL regime near 1°, with typical CDM models having a small 
amount of linear bias {b\ ~ 1.2 ± 0.1). The thin horizontal line is taken from Figure |3|. 
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itself an unbiased estimator (Hui & Gaztafiaga 1998), and 3) nonlinear effects not accounted 
for in our models (Boucher, Schaeffer, & Davis 1991; Jain & Bertschinger 1994; Scoccimarro 
et al. 1998; Scoccimarro & Frieman 1998; Munshi & Melott 1998). At large scales we 
expect discrepancies to arise from the breakdown of the small-angle approximation and the 
large scatter in the data induced by the smaller number of independently sampled cells. If 
we therefore exclude the data points at both the lower and upper ends of the range of 0, 
and restrict the comparison to QL scales of 6 < 1°, where our predictions are expected to 
be valid, we find reasonable agreement with plausible CDM models having a small amount 
of linear bias. Of course, this conclusion is merely qualitative, as these data are clearly 
unable to place tight constraints on the models. The 2-er error on b\ in both models is 
roughly ±0.1, depending on exactly which points are excluded, in agreement with previous 
determinations of b\ ~ 1.0 ±0.2 for the APM Survey (Baugh & Efstathiou 1993; Gaztafiaga 
1994; Bernardeau 1995). That the data for optically selected galaxies are roughly fit in 
the QL regime by the predicted values of S3 in CDM models, requiring little or no bias 
(Gaztafiaga 1994; Bernardeau 1995; Cappi & Maurogordato 1995; Roche & Eales 1998; 
Magliocchetti et al. 1998), lends some support to the scenario of gravitational evolution 
ensuing from Gaussian initial conditions, but better data are clearly needed. Note that 
the thin solid line, derived from equation (21) with the APM SF, no bias, and a constant 
power-law index of n — — 1.2, yields a poorer fit in the QL range. 



3. ANGULAR THREE-POINT CORRELATION FUNCTION 

The normalized angular skewness, while containing some information about the overall 
dependence of clustering strength with scale, discards detailed and valuable information 
about the configuration dependence which is contained in the full angular 3PCF. It is 
known that the shape dependence of the spatial bispectrum and 3PCF, in leading order 
PT, can be used to distinguish between the effects of gravitational clustering and bias, as 
well as to obtain independent constraints on the linear- and nonlinear-bias parameters and 
on cosmological parameters (Jing & Zhang 1989; Fry 1994; Fry 1996; Jing & Bonier 1997; 
Matarrese, Verde, & Heavens 1997; BK99). It stands to reason that the full angular 3PCF 
should also break these degeneracies and yield similar constraints. With this in mind, many 
authors have investigated the normalized angular 3PCF, 

„fa a a \ _ Z tglg^ag, #3l) /o 9 \ 

91^12,^23,^31; — 7^-^ 77 \ , 777— \ 77 \ , TTf^ 77 T, 

w(6 12 )tv{623) + w{9i 2 )w{e u ) + 07(6*23) 07(6*31) 
where the angular distances 6V,- form a triangle on the sky. 



This quantity, which like S3 is independent of the amplitude of the power spectrum, 
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was introduced based on empirical arguments for the so-called hierarchical model (Peebles 
1975; Peebles & Groth 1975; Peebles 1980), and motivated by data which appeared to 
indicate that q was simply a constant. While equation ( |3"2] ) may be taken as a definition, 
the result for q in QL PT is, in general, not merely a constant; we demonstrate that 
the apparent constancy of q, as measured from past data, may simply be an artifact of 
coarse averaging over configuration shapes (Peebles & Groth 1975; Fry & Seldner 1982), 
combined with the effects of the shallow SFs which have characterized past surveys. There 
are, however, other problems associated with the definition of the normalized angular 
3PCF in equation (|32|). In particular, q is apt to exhibit rapid variation and divergence 
where the denominator in equation ( p2|) happens to acquire values of or near zero.[] This 
behavior is artificial in the sense that it does not reflect the true behavior of the 3PCF 
itself. BK99 demonstrate the analogous problems associated with the definition of the 
normalized spatial 3PCF, Q, and speculate that practical limitations in measuring this 
behavior may, in part, be responsible for observed discrepancies between QL PT predictions 
and those of iV-body simulations. The skewness is less susceptible to this problem since 
it is scaled by the square of the variance, a positive-definite quantity. In places where q 
exhibits this rapidly- varying or divergent behavior, we will instead consider qy, defined 
as qy = Z(6'i2, ^23, 03i)/(aio) 4 , where a\° = ^Jzuio(0) is the rms angular fluctuation on a 
smoothing scale of 1°, obtained for each model using equation (pTTj) . The chosen smoothing 
scale is, of course, arbitrary; we adopt a value of 1° as a convenient choice for the sake 
of graphical representation. Like q, qy is independent of the overall normalization of the 
power spectrum, but it is not susceptible to the above mentioned problem associated with 
the definition of q. Observational measurements of both q and S3, however, will necessarily 
be ratios of estimators, and thus subject to associated errors discussed by Hui & Gaztanaga 
(1998). 

Taking equation (^) with x = k8 and omitting the window function, since we presume 
that the full CFs are now measured by direct counting, rather than by smoothed counts in 
cells, we obtain the direct-counting angular 2PCF, 

The full angular 3PCF for projected triangles with sides of angular measures #12, #23, and 



7 We note that this problem is avoided entirely if one considers instead the normalized angular bispectrum 
in Fourier space. 
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6*3i can be evaluated using equation (|T^ 
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where 0j is the angle between re, and 0^3, -0 is the angle between K\ and k 2 , and similar 
conventions are used in each permutation term. Taking = ftj^, and 0i 2 to be the angle 
between #i 3 and #23, so that 

cos ip = cos 0i2 (cos 0i cos 02 + sin 0i sin 02 ) + sin 0i2 (cos 0i sin 2 — sin 0i cos 02) , (35) 

and performing the 0i and 02 integrations, we have 
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cos 2 012 J 2 (xi) J 2 (x 2 ) + -W{xi)W{x 2 ) 



-~J 2 ( Xl )W(x 2 ) - ~ J 2 {x 2 )W{x x ) 



(eye). 
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For evolving bias we simply substitute 



1 



C 2 (w), 



1 — 



C 3 H 



(37) 



inside the brackets, and replace &i by 1/C2(w). 



Equation (^) does not contain terms analogous to the higher-order terms in the sum 
in equation fl2~0|) [which gave rise to the A(n) term in equation (21)], since evaluation of 
the 3PCF by direct counting requires no smoothing. Using our model for unsmoothed bias, 
such terms would be induced in the calculation of the 3PCF using counts in cells, arising 
from factors of the form W(xQ/9) in the appropriate integrands. One must be careful to 
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realize that if the unsmoothed-bias prescription is correct, the nonlinear-bias parameter 
measured from the direct-counting angular 3PCF in equation fl3"6|) will be different from the 
smoothed value which could be inferred from S3 or from the counts-in-cells 3PCF, though 
the two can be related, e.g., via the A(n) term in equation (fH|). 



3.1. Effects of the Power Spectrum 

For a scale-free power spectrum, equation (|33|) yields the familiar result, w{9) oc 0~( n+2 \ 
For realistic CDM power spectra, there is a further dependence on 9 which arises from the 
scale-dependence of the transfer function. Likewise, had we ignored the shape-dependent 
terms in the bispectrum in equation (fToT), we would find for power-law spectra the 
hierarchical approximation Z{9\2, 023, #31) oc w(0i2)^(023) + vo{9i2)w{9zi) + w (6*23)^(^31 )■ 
These terms, however, introduce into equation (|36"D an additional, non-trivial dependence 
on the shape of the triangular configuration, even for the case of scale-free power spectra. 
Realistic models must account for this variation, as well as the scale dependence induced 
by the transfer function. 

In order to interpret the results of the various CDM models, it will be instructive to 
first investigate the results obtained for Z for scale-free power spectra. In this case, we have 



£(012,023,031) = Rs(n) 



10/3 / 013 023 \ 

w(9 13 )za(9 23 ) + (eye), (38) 



+ 4 (2 - 2(3 + (3 2 cos 2 12 



7 (2-/3) 2 

where w(9) oc _/3 and (3 = n + 2. Figure || displays the results for q/R^{n) for 
n = { — 1,-1.5,-2} [recall R-^n) is of order unity, equal to {1.26, 1.19, 1.15} and {1.21, 
1.18, 1.17} for these values of n, for the APM and high-z SFs, respectively]. Since q in this 
case is independent of overall scale, we specify the triangle configuration by fixing # 13 to 
have length 1, #23 to have a relative length u > 1, and examine the variation of q with 0i2, 
for u = 1,3, and the "limit" of u — > 00. The Figure reveals shape dependence on both u and 
0, the extent of which can (particularly in the case of the dependence) depend strongly 
on the index n. Specifically, we note that the amplitude of the shape variation tends to 
increase with increasing n. An "average" value for q can be obtained by taking (cos0) = 
and (cos 2 </>) = 1/2 in equation (j38[), yielding q = yi?3(n); while this is not generally 
valid, since the three angles appearing in the cyclic terms cannot be varied independently, 
this result is exactly true for all configurations when (3 = (the no-smoothing case), as 
demonstrated by the thin, solid, horizontal line. 
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Fig. 8. — Results for q(9i S , 6* 23 , <j)vi)/Rs{n) for power-law spectra with indices of n — — 1, 
(thick lines), —1.5 (thin lines), and —2 (thin horizontal line). Since these results are 
independent of overall scale, we have fixed #13 to have length 1, and vary #23 in the manner 
shown. Note that the dependence of q on triangle shape varies strongly with n. 
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Fig. 9. — The left panels show results for q, assuming the APM SF, using the {#13, 623, ^12} 
parameter set to define triangles. The right panels show the results for qy for the same 
choice of configurations (assuming 12 > 60°) using the {r,u,v} parameter set. In both 
cases, the unbiased SCDM (solid lines), TCDM (long-dashed lines), ACDM (short-dashed 
lines), and OCDM (dotted lines) models are shown. Note that for both SFs, models with 
different values of T yield distinct predictions, and that the precise configuration dependence 
of the normalized 3PCF in each different model can vary appreciably with the SF. 
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In the left panels of Figure ||| we illustrate, for the APM SF, the behavior of 
9(^13,^23,012) for several configurations falling in the quasilinear regime, for an unbiased 
tracer mass in the SCDM, TCDM, ACDM and OCDM models. In this case, the results 
for the latter two models differ only by about the thickness of the curves, and can barely 
be distinguished. Other models, however, yield quite distinct predictions and might be 
differentiated on the basis of angular 3PCF measurements from large surveys. In particular, 
the SCDM and TCDM can be clearly distinguished from each other (BKJ show that this 
is not the case when considering the spatial 3PCF at z — 0), and from the low-r models, 
which show comparatively little variation with 12 . Generally, we now find a dependence of 
q on scale (compare, e.g., the top and bottom panels) arising from the shape of the CDM 
transfer function. The observed shape dependence on 6*23 (compare, e.g, the top and middle 
panels) and i2 of these models can be understood in terms of the results for scale-free 
power spectra: over the scales considered here, the SCDM model has the largest effective 
power-spectrum index, while the low-r models have the smallest. The results for the APM 
SF shown in Figure [| roughly agree with the predictions and iV-body results of Frieman 
and Gaztanaga (1999). We thank them for pointing out an error in an earlier calculation of 
ours. 

For a shallow SF such as that describing the APM Galaxy Survey, three points 
describing an elongated triangle in projection are more likely to correspond to a 
configuration which is elongated in real space. Thus, our results for q using the APM SF 
are similar to corresponding results for the normalized spatial 3PCF, Q (Fry 1984; Jing & 
Borner 1997), which predict that clustering in the mildly nonlinear regime favors elongated 
structures. For low-r models, however, the APM SF happens to yield predictions which 
do not vary strongly with either the shape or size of the triangular configuration. Since it 
is known that measurements of the power spectrum from surveys such as the APM favor 
models with low values of T (Efstathiou, Sutherland, & Maddox 1990; Baugh & Efstathiou 
1993), it is not surprising that measurements of the angular 3PCF in shallow surveys have 
indicated it to be configuration independent, particularly given the present observational 
errors (Peebles and Groth 1975; Groth & Peebles 1977; Jing & Zhang 1989; Toth, Hollosi, & 
Szalay 1989; Jing, Mo, & Borner 1991; Borgani, Jing, & Plionis 1991; Jing, & Borner 1998). 
While this has no doubt has lent support to the empirically-motivated hierarchical model 
which assumes Z to be, in fact, shape independent (i.e., q =constant), we demonstrate 
that the angular 3PCF, in general, can depend more significantly on the configuration and 
analysis of deeper angular surveys must take this variation into account. 

The results for high- z SF for the same configurations (assuming 0i 2 > 60° are shown in 
the right panels of Figure |9|. We now plot qv, since for this choice of SF the combination of 
2PCFs appearing in the denominator of q can acquire values of or near zero in our models. 



Furthermore, we now take the defining parameters for triangles projected on the sky to be 



r = 9 



12, 




V = 




(39) 



where 9 12 < #23 < #31, so that u > 1 and < v < 1 (Peebles & Groth 1975). For this choice 
of parameters, r effectively fixes the overall size of the triangle, while u and v determine the 
exact size and shape. Figure 6 in BK99 provides a graphical illustration of the variation of 
triangle geometry with r, u, and v. This parameterization is convenient since each choice 
of {r, u, v} involves a minimum scale set by r. The {#13, #23, 4>u} parameter set, with 
9 13 = r and #23 — ur is equivalent for 12 ranging down to 60°, but for less than 60° this 
parameter set can sample scales (given by 9 12) all the way down to zero. Since we have 
chosen to normalize Z by the square of the variance at a fixed angular scale, the {#13, #23, 
012 } set can yield dramatically different amplitudes for q v if 0° < </> 12 < 60° as opposed to 
those for 60° < 0i 2 < 180°, and is therefore less tenable. 

Even once their different normalizations have been accounted for, we find that the 
high- z curves have significantly lower amplitudes than their APM-SF counterparts in the 
left panels. This is simply because projected triangles in a deeper survey can correspond 
to much larger physical triangles, for which the clustering amplitude will be significantly 
smaller. Since qv is normalized by the square of the variance at a fixed scale, these plots 
now reveal clearly the expected decrease in clustering strength with increasing scale; this 
behavior, however, is not described simply by the Z oc 9~ 2n+A prediction of the hierarchical 
model. Unlike the APM case, we now find that all models show a significant variation with 
shape, as well as scale, with the low-r models in some cases now showing the strongest 
shape dependence. The SCDM and TCDM models appear less differentiated than in the 
APM case, but we find that the ACDM and OCDM models can now differ somewhat 
more substantially, particularly on larger scales, since deeper surveys are more sensitive 
to cosmological effects. Most importantly, however, we still find overall that the high-r 
and low-r models can be fairly well resolved on 2° scales, and the distinction increases 
for larger scales. Naturally, a thorough comparison with measurements spanning the full 
configuration space in a more complete fashion would allow for better discrimination among 



models. 



3.2. 



Effects of Bias and its Evolution 



The results for the different bias scenarios investigated in §2.2 are demonstrated in 
Figure [II], which shows qv, assuming r = 1° and u = 1, for the ACDM (thick lines) and 
OCDM (thin lines) models, for the high- z SF. Though this is only one slice of the parameter 
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Fig. 10. — Plot of qv(v), assuming fixed r = 1° and u = 1, for the five bias scenarios from 
§2.2, each employing the ACDM (thick lines) and OCDM (thin lines) models, and using the 
high- z SF. Note the respectively similar influences of the biasing models as were seen with 
s 3 (Figures || and ||). 
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space we are exploring, we use it to illustrate the similar dependences of the full angular 
3PCF on the bias scheme as were seen for S3. In particular, we see a general flattening and 
reduction of qy with increasing linear bias, and a relative increase with increasing nonlinear 
bias, as compared with the corresponding linear scheme. Note the slightly different 
predictions between these two models, arising from their different expansion histories; these 
differences would be greater for lower values of fl . As with s 3 (see Figures ^| and H), it is 
clear that measurements of qy should be able distinguish between the effects of a significant 
bias and those arising simply from unbiased gravitational evolution. 

Consideration of the full configuration dependence of the angular 3PCF, however, 
allows for much better discrimination between different cosmological and biasing schemes 
than could be obtained from the skewness alone. It is known that measurements of the 
bispectrum and 3PCF can be used to obtain independent constraints on the constant- valued 
linear- and nonlinear-bias parameters (Fry & Gaztanaga 1993; Fry 1994, 1996; Jing 1997; 
Matarrese, Verde, & Heavens 1997, BK99). This can be inferred here from the different 
manners in which b\ and 62 appear in equation (|3~6D ; the different weightings of these two 
parameters is sampled for every triangular configuration, producing a graph such as Figure 



10| for every pair of values (r,u) [as opposed to s 3 , which yields only one constraint on 
the combination of b% (t) and 62 (^)]- As another example of the added information which 
can be gleaned from the full geometric dependence of the full 3PCF, we plot in Figure 
|TT| q v (r,u,v), assuming the high- z SF, for two of our bias scenarios: the SCDM constant, 
linear-bias model (with b\ = 2 and b 2 = 0), and the SCDM evolving linear- bias model (with 
&i„ = 5 and 62, = 0). These two cases yielded nearly identical results for S3 in Figure 
^| (where evolving, linear bias effectively approximated a smaller, constant, linear bias), 
particularly on scales near 1°. Here, we see that by considering the variation of qy with 
triangle geometry, the degeneracy between constant and evolving bias can be alleviated, to 
a degree, because the angular 3PCF probes a broad range of redshifts and allows us to see 
the effects of evolution in projection. BK99 show that these evolving and non-evolving bias 
models might also be distinguished using measurements of ( as a function of redshift in a 
deep (z <: 1) survey, or, to a lesser extent, by large-scale (R ^ 20 Mpc) measurements of ( 
at z = 0. In practice, however, measuring the angular 3PCF over scales of ~ 1° in a deep 
optical survey, or a radio survey such as the VLA FIRST Survey, may be more feasible and 
would likely yield stronger statistical constraints. 

The parameter values we have investigated here are only intended as illustrations. 
The full statistical power which could be gained by exploring both the configuration and 
evolutionary dependences of Z in a complete fashion could be useful in distinguishing 
between the signatures of gravitational evolution versus bias, independently measuring the 
linear- and nonlinear-bias parameters, and constraining the shape of the power spectrum. 
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Fig. 11. — Plot of qv(r,u,v), assuming the high- z SF, for the SCDM model with constant, 
linear bias (dotted lines) and the SCDM model with evolving linear bias (short-dashed lines) 
from Figure |6|. While these models yielded nearly degenerate predictions for S3 in Figure 
H, they are better distinguished here by virtue of the different configuration dependences of 
their full angular 3PCFs. 
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With better data, as might be obtained from emerging surveys, one could further determine 
the linear evolution of the power spectrum and of the bias, which in turn provide direct 
constraints on Q , Q\, and possibly the epoch of galaxy formation. However, it is already 
clear from Figures |9|-|TT| that the behavior of the angular 3PCF is not simply that of the 
hierarchical model which assumes q to be a shape-independent constant, and further that, 
in addition to the full geometric dependence, one must consider the detailed dependences on 
the forms and evolutionary behaviors of the power spectrum, the biasing mechanism, and 
the selection function, as these can have a substantial impact on the results, particularly for 
high-redshift surveys. 



4. CONCLUSION 

We have calculated the leading-order results for the normalized angular skewness 
and 3PCF, assuming Gaussian initial conditions, for an arbitrary, biased tracer-mass 
distribution in flat and open universes. We have considered two selection functions, 
intended to represent low-redshift (z ~ 0.1) and high-redshift (z ~ 1.0) surveys. Our 
derivations incorporate such features as the scale dependence and linear evolution of the 
CDM power spectrum and the presence of a possibly evolving linear or nonlinear bias. 
Our results for S3 for the case of an unbiased tracer mass are in agreement with previous 
work (Fry & Gaztahaga 1993; Bernardeau 1995) and we confirm that PT predictions 
for S3 in typical CDM models, with a small amount of linear biasing (bi ~ 1.15 ± 0.1), 
agree in the QL regime with data from the APM Galaxy Survey. Similar to BK99, 
we extend the result for S3 in the case of nonlinear bias, to include a scale-dependent, 
leading-order correction which becomes significant for positive effective spectral indices, 
corresponding to scales > 5° for CDM models. This correction term can be used to 
relate the unsmoothed nonlinear-bias parameter defined in the limit of continuous fields, 
with the traditionally-used, scale-dependent, nonlinear-bias parameter defined for smoothed 
density fields. For a given CDM model, this large-scale behavior of S3 might alone allow 
a more accurate determination of the linear- and nonlinear-bias parameters (or at least 
discriminate between the smoothed and unsmoothed bias formulations), if other factors, 
such as sampling variance and the breakdown of the small-angle approximation, can be 
addressed. 

We also derive predictions for the normalized angular 3PCF, considering qy rather 
than q in cases where the latter is ill-behaved. We show that Z, as predicted by QL 
PT, is not given simply by the form of the "empirical" hierarchical model, i.e., with q 
simply a constant of order unity, but has a more complex dependence on scale, shape, 
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and expansion history. While for shallow SFs the weak shape dependence of q, combined 
with measurements coarsely averaged over configuration shapes, can yield results consistent 
with this hierarchical model, ignoring the detailed geometric and evolutionary variations 
of the angular 3PCF in general can destroy much valuable information. In several cases, 
we illustrate explicitly the importance of the various evolutionary dependences which 
factor into our calculations of angular statistics, such as the linear evolution of the power 
spectrum, the evolution of the bias parameters, and the f2 dependence of the SFs. The 
proper treatment of these detailed dependences will be important in interpreting data from 
the emerging generation of deep, high-precision surveys, particularly when considering the 
high-redshift universe. 

We find that biasing plays a major role in predictions for s 3 and Z. In particular, the 
scale dependence of S3 and the configuration dependence of q (or qy) bear characteristic 
imprints of bias, such as a relative flattening and decrease with increasing b\ and a relative 
increase with increasing 62- Compared with the dependence on the biasing scenario, the 
dependence of S3 on the cosmological model is fairly weak. This is not the case, however, for 
q, which can vary appreciably with the values of T and n. Though S3 preserves information 
about the overall scale dependence of the 3PCF, it provides only one constraint on the 
combination of the cosmological parameters, bi(t), and 62 (t), leaving various possible 
degeneracies between these quantities, including the degeneracy between constant- and 
evolving-bias models. These can be partially alleviated by considering the dependence 
of the normalized angular 3PCF on the projected-triangle geometry. In particular, the 
variation of q can be used to distinguish between the effects of gravitational evolution and 
bias, place constraints on the value of T and the shape of the power spectrum, measure the 
bias terms b\ and 62 and their time evolution, and perhaps constrain Q directly. Comparing 
predictions for S3 and q with several data sets characterizing different tracer populations 
could allow multiple, complementary constraints on these important quantities. 

While many of the results obtained here for S3 and Z are qualitatively similar to the 
corresponding results obtained by BK99 for their spatial counterparts, S3 and (, the use of 
angular statistics offers many practical advantages, such as greater statistical power and the 
absence of redshift distortions, which warrant their separate consideration. Our derivations 
have also relied on specific assumptions, such as the chosen bias-evolution model and 
the validity of both the small-angle approximation and the assumption of quasilinearity. 
Based on analogous results from three-dimensional calculations, we expect that qualitative 
features such as the upturn in s 3 at large smoothing angles in the case of nonlinear bias, 
the reduced configuration (scale) dependence of q (S3) with increasing linear bias of any 
form, and the sensitivity of q to the shape of the power spectrum will be retained in more 
general calculations. 
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The utility of the predictions for S3 and q will of course rely on the precision with 
which these statistics can be measured. While current data, of the type shown in Figure [7|, 
is clearly insufficient to obtain strong constraints on the shape of the power spectrum or on 
the bias parameters, it is hoped that future surveys, such as the SDSS, will greatly improve 
upon this situation. In particular, if we optimistically assume the relevant uncertainties 
to be predominantly statistical, we can infer that the errors in measurements of S3 from 
the SDSS will be roughly a factor of ten smaller than those in Figure [7]. Of course, the 
comparison of predictions with two-dimensional survey data will require the additional 
consideration of systematic effects. Among these are finite-volume and boundary effects, 
which can have a significant impact on degree scales and above (Baugh, Gaztanaga, & 
Efstathiou 1995; Szapudi & Colombi 1996; Colombi, Szapudi, & Szalay 1998; Gaztanaga 
& Bernardeau 1998; Munshi & Melott 1998), estimation biases (Hui & Gaztanaga 1998; 
Kerscher 1998), Poisson noise (Peebles 1980; Gaztanaga 1994), and sampling variance (Seto 
& Yokoyama 1998). Gaztanaga & Bernardeau (1998) also find that QL PT predictions for 
angular statistics may underestimate projection effects, particularly for APM-type surveys. 

In addition, the above calculations can be extended in several ways. For example, one 
might consider stochastic-biasing models which may yield different results, particularly on 
smaller scales (Fry & Gaztanaga 1993; Catelan et al. 1998; Catelan, Matarrese, & Porciani 
1998; Taruya, Koyama & Soda 1998; Blanton et al. 1998; Colin et al. 1998; Taruya & 
Soda 1998; Narayanan, Berlind, & Weinberg 1998). Several authors have also investigated 
the impact of higher-order nonlinear terms in PT calculations (Jain & Bertschinger 1994; 
Scoccimarro & Frieman 1996a,b; Scoccimarro et al. 1998; Heavens, Matarrese, & Verde 
1998; Scoccimarro & Frieman 1998), also expected to become important on smaller scales, 
where w > 1. Similar calculations can also be performed for higher-order angular moments 
and n-point CFs, such as the kurtosis and four-point correlation function, and so on 
(Suto & Matsubara 1994; Bernardeau 1994; Lokas et al. 1995; Gaztanaga & Bernardeau 
1998), which, in principle, provide additional, independent constraints. The formalism 
can also be applied to structure-formation models with non-Gaussian initial conditions, 
such as topological-defect or isocurvature models, to investigate how these might be 
distinguished from inflationary models (Luo & Schramm 1993; Jaffe 1994; Fry & Scherrer 
1994; Chodorowski & Bouchet 1996; White 1998). 

We are indebted to J. Frieman and E. Gaztanaga for helpful discussions and suggestions, 
and in particular for their help in clarifying an earlier problem in the calculation. We also 
wish to thank R. Scoccimarro, as well as Licia Verde, for numerous helpful comments and 
discussions. This work was supported at Columbia by D.O.E. contract DEFG02-92-ER 
40699, NASA NAG5-3091, NSF AST94-19906, and the Alfred P. Sloan Foundation, and at 
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Berkeley by NAG5-6552. 

A. PROJECTED BISPECTRUM 

In this Appendix, we derive a simple Fourier-space analogue of Limber's equation for 
the bispectrum of the two-dimensional projection of a three-dimensional scalar field. Our 
derivation is similar to that of Kaiser (1992), who calculated the projected (two-point) 
power spectrum. We first consider the case of an Einstein-de Sitter universe and at the end 
we generalize the results to open or closed universes. Jaffe & Kamionkowski (1998) discuss 
a different application of this formalism to isotropic vector fields. 

Consider a three-dimensional field, 5(x,w), at a conformal lookback time w, and 
its Fourier transform, 5(k,w). We will take 5 to have zero mean (or subtract the mean 
explicitly). To project onto the two-dimensional sky, we apply the "selection function," 
q(w), and denote the resulting field by 

POD 

p{6) = / dw q(w)5(w9i,w9 2 ,w,w), (Al) 
Jo 

where 6 is a two-dimensional vector on the sky pointed at by the three-dimensional unit 
vector 9. We consider only small patches, and thus can perform a two-dimensional Fourier 
transform on essentially flat areas of sky, giving p(k,), where k, denotes two-dimensional 
Fourier-space vectors. As in Kaiser (1992), the contribution from a thin shell of width Aw 
centered at wo to p(k) is 

~ Awq(w ) (dk^- fKx k 2 \ . (k 3 Aw\ 

Ap{K) — 2 / V S — >— M,w)3q , A 2 

Wq J 2n \wo w J \ 2 J 

where the jo{x) = sinx/x factor arises from the integral of exp(— iks • w) over the shell. 

We define the three-dimensional bispectrum, the Fourier-space analogue of the 
three-point function, by 

(S(k)5(k')5(k")) = (2n) 3 5 D (k + k' + k")B(k, k', k")- (A3) 

the two-dimensional analogue is 

(P(k)p(k')p(k")) = (2tt) 2 5 d (k + k' + k")B p (k, k', k"), (A4) 



where the dimensionality of the Dirac delta function is obvious from its arguments. 
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We can form the contribution to the two-dimensional bispectrum from the thin shell, 

'k 3 Aw\ . fk' 3 Aw\ . {kg Aw* 



Wq J 2n 2tt 



x io^^-)io^^-Jio^^-JS(r,r',r"), (A5) 



where 



\ — r- + *§■ (A6) 

^ w o 

Now, we make the same simplifying assumptions as in the original Kaiser derivation: the 
shells have width Aw/w <C 1 but are thick compared to wavelengths contributing to S of 
interest, so Aw/w 3> 1/k. Since the jo factors have width Ak 3 ~ I /Aw, we can pull the B 
out of the integral and set r = r| fe3=0 = k/w ; only modes perpendicular to the line of sight 
contribute in the small-angle approximation. Then using the identity, 

OO POO 

2 



3o(x) j (y) j (x + y)dxdy = it , (A7) 
we can write the contribution from the shell as 



oo J— oo 



A P (k)A P (k')A P (k")) = (2n) 2 5 D (K + k' + K '') q ^ w °\ Aw B ( -, -, — ) . (A 



Wq \W W W 

The bispectrum is a "cumulant" ; that is, contributions to the total bispectrum simply add. 
Thus, we can sum up the contributions from different shells, and convert the sum to an 
integral: 

(p(k)p(k')p(k")) = (27t) 2 5(k + k,' + k") f dw ^^B ( -, -, — ; w] ; (A9) 

J w* \w w w J 

with the full definition of the bispectrum, this just gives the desired formula, 

B p (k, k', k") = fdw ^B (-, -,-;w) (AlO) 
J wr \w ww J 

(allowing the three-dimensional bispectrum to vary with lookback time, w). Note that the 
integrand in the formula is weighted by q 3 /w A ; this is in contrast to the (two-point) power 
spectrum version of the calculation which weights by q 2 /w 2 [see equation (|^)]. Hence, the 
bispectrum is weighted considerably more heavily to late times (small w). 

To generalize to an open or closed Universe, we must replace the distances w in the 
derivation above with the angular-diameter distance. Thus in an open or closed Universe 
one replaces 

1 clqH \J\1 — Q — Q\\ 

> / (AH) 

w S(a H wJ\l-Q -Q A \) 
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in the first three arguments of B (k/w, k'/w, k"/w; w) and in the w 4 in the denominator 
in equation ( |A10| ), where S(x) = sinhx in an open Universe and S(x) = sinx in a closed 
Universe. 
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